^^5851 0 


SENSITIVITY  ANALYSIS  IN 


ULTITERMINAL  FLOW  NETWORKS 


by 

Francr*co  K.  Rack) 


OPERATIONS  RESEARCH  CENTER 

COLLEGE  OF  ENGINEERING 


UNIVERSITY  OF  C  A  L  I  F  0  R  N  I A  -  B  E  R  K  E  L  E  Y 


P 


SENSITIVITY  ANALYSIS  IN  MULTITERMINAL  FLOW  NETWORKS 


by 


Francisco  K.  Rado 
Operations  Research  Center 
University  of  California,  Berkeley 


June  1967  ORC  67-22 

This  research  has  been  partially  supported  by  the  U.S.  Army  Research 
Of f ice-Durham,  Contract  DA-31-124-ARO-D-331 ,  and  the  National  Science 
Foundation  under  Grant  GK-1684  with  the  UnLversity  of  California. 
Reproduction  in  whole  or  in  part  is  permitted  for  any  purpose  of  the 
United  States  Government. 


Algorithms  are  given  for  determining  the  influence  of  varying 
the  capacity  of  an  arc  or  of  several  arcs  with  linearly  related 
capacities,  on  the  maximal  flows  between  every  pair  of  nodes  in 
an  undirected  network.  This  work  builds  upon  the  fundamental 
work  of  Gomory  and  Hu,  and  is  closely  related  to  the  approach 
of  Elmaghraby,  who  also  gives  algorithms  for  solving  these 
problems. 


SECTION  X 


Description  of  the  Problem 

0 

Given  a  network  G  In  which  all  arcs  have  finite  capacities,  a 
standard  computational  procedure  due  to  Ford  and  Fulkerson  [1]  yields 
the  maximal  flow  from  one  node  called  the  source  to  another  called  the 
sink. 

The  problem  of  determining  the  maximal  flows  between  every  pair  of 
nodes  in  an  undirected  network  Is  called  the  maximal  flow  multi-terminal 
problem.  Although  there  may  be  many  pairs  of  nodes,  there  are  at  moat 
n-1  different  maximal  multi- terminal  flows,  where  n  la  the  number  of 
nodes  In  G  . 

The  first  problem  solved  In  the  present  work  Is  that  of  determining 
how  varying  of  the  capacity  of  an  arc  In  a  given  range  affects  the  n-1 
values  of  the  maximal  flows.  The  second  one  generalizes  the  first  by 
allowing  the  capacities  of  any  subset  of  arcs  of  the  network  to  depend 
linearly  on  a  parameter  constrained  to  a  specified  range,  and  then  asks 
for  the  relation  between  the  parameter  values  and  the  maximal  flow  values. 

Briefly,  the  first  problem  is  the  sensitivity  analysis  of  multi- 
terminal  flow  networks  to  changes  in  the  capacity  of  one  arc.  The  second 
extends  this  analysis  to  several  arcs  for  a  special  kind  of  network;  namely, 
one  In  which  the  capacities  of  the  arcs  are  linearly  related.  Our  chief 
result  Is  the  derivation  of  two  algorithms  that  solve  these  problems. 
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Introductory  Section 

Because  the  fundamental  results  concerning  maximal  flows  In  networks 
are  well  known  (see,  e.g.,  [1]),  we  will  review  them  here  only  to  the 
extent  necessary  to  establish  the  notation  and  definitions  useful  to  this 
paper.  Let  G  ■  G(N,A)  denote  a  network  where  N  Is  the  set  of  nodes 
and  A  the  set  of  arcs  of  G  .  The  arc  a^  joining  node  1,  to  node  £ 
has  capacity  If  the  maximum  amount  of  flow  that  this  arc  can  carry 

from  node  1,  to  node  J.  Is  C.  . 

» 

If  the  arc  can  carry  a  flow  of  units  either  from  _1  to  j  or 

from  j,  to  1  it  is  said  that  the  arc  is  undirected.  If  every  arc  of  A 
has  that  property,  we  speak  of  an  undirected  network. 


Maximum  flow  -  Minimum  out  Theorem t  Corollaries  and  algorithms 
related  to  it: 


A  cut  C  in  G(N,A)  separating  nodes  £  and  t  Is  a  set  of  arcs 
(a^j/leX,  JeX)  ■  (X,X)  where  £  e  X  and  £  e  X  .  It  also  partitions  the 
nodes  In  two  subsets:  X  and  X  .  The  capacity  C(X,X)  of  the  cut 
(X,X)  is  the  sum  of  the  capacities  of  the  arcs  of  (X,X)  . 

The  maximal  flow  -  minimum  cut  theorem  states: 

For  any  network  che  maximal  flow  £  from  £  to  £  is  equal  to  the 
minimal  cut  capacity  of  all  cuts  separating  £  and  t  . 

A  path  from  node  X^  to  node  XR  la  a  sequence  of  arcs 
{(X^,  X^)},  k-1,  •••.  n-1  ,  where  X^  4  Xp  if  kjlp  . 

Some  corollaries  of  the  theorem,  for  undirected  networks,  are: 

A  cut  (X,X)  is  minimal  if  and  only  if  every  maximal  flow  f_  saturates 
all  arcs  of  (X,X)  with  flows  going  from  X  to  X  .  A  flow  £  Is  maximal 
If  and  only  if  there  is  no  flow  augmenting  path  with  respect  to  £  . 


The  algorithms  for  this  problem  are  systematic  searches  of  paths  from 
source  to  sink,  such  that  through  each  path  we  can  carry  a  new  flow  given 
by  the  minimum  of  the  differences  of  the  capacities  of  the  arcs  of  the  path 
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and  the  flows,  In  the  direction  of  the  path,  already  carried  by  them. 

Multi-terminal  network  flews: 

Sometimes  the  knowledge  of  the  maximal  flow  between  several  or  all 
pairs  of  nodes  can  be  very  useful.  By  naming  any  2  nodes  ±  and  J. 
source  and  sink,  we  can  find  the  maximal  flow  f^  between  ^  and  by  the 
standard  algorithms.  Applying  this  procedure  j  (n-1)  times,  (the  number 
of  conAlnations  of  two  nodes  In  n)  yields  a  maximal  flow  between  every 
pair  of  nodes.  (The  maximal  flow  limited  by  capacity  between  JL,  and  J[ 
is  the  same  as  the  maximal  flow  between  J[  and  .1  for  an  undirected 
network.) 

Gomory  and  Hu  [3]  have  shown  that  there  are  not  more  than  n-1  different 
maximal  flows  and  they  .an  be  found  by  solving  only  n-1  maximum  flow 
problems  (some  of  which  are  determined  in  smaller  networks  than  the 
original) . 

To  describe  the  Gomory  and  Hu  approach,  we  Introduce  the  following 
additional  definition: 

A  connected  set  of  arcs  Is  one  in  which  there  is  a  path  between  every 
two  nodes . 

A  ayole  is  a  path  in  which  the  first  and  last  node  coincide. 

A  tree  Is  a  connected  set  of  arcs  with  no  cycles. 

A  spanning  tree  of  a  network  Is  a  tree  that  Includes  all  the  nodes  of 
the  network. 

A  maximal  spanning  tree  is  one  In  which  the  sum  of  the  capacities  of 
the  arcs  Is  maximal. 
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He  now  show  that  thcra  are  at  most  n-1  different  values  of  f^  ,  (l.jeN) 

Theorem  (2):  A  necessary  and  sufficient  condition  for  all  f^  , 

*Jk'  *lk  to  be  the  max^num  flows  between  nodes  ,1  and  ±  ;  ±  and  k  ;  and 
_1  and  k  respectively  la  that 

(1)  «uim„  (fir  fJk)  . 

To  short  necessity:  Suppose  f^k  <  Min  (f^,  fjk)  •  then  by  the  maximal  flow- 
minimal  cut  theorem  there  Is  a  cut  (X,  X)  with  node  1  in  X  and  k  In 
X  ,  such  that  the  capacity  of  the  cut  la  equal  to  f^k  .  Now  node  J.  belongs 
either  to  X  or  to  X  . 

If  it  la  in  X  ,  fJk  >  C(X,  X)  . 

It  it  la  in  X  ,  f  >  C(X,  X)  . 

In  both  cases  a  contradiction  Is  obtained. 

(2)  By  induction  fip  >  Min  (f^,  fjk,  ffcl,  ...»  f  )  . 

To  short  sufficiency:  Let  H  a  "modified"  network  composed  of  the  same  nodes 
G  but  with  arcs  connecting  every  pair  of  nodes  1,  j  and  with  capacities 
equal  to  the  f^'s  of  G  . 

Now  we  construct  the  maximal  spanning  tree  of  H  .  We  shall  refer  to 
the  arcs  of  this  tree  as  "links."  Any  f^  of  H  that  does  not  appear  as  a 
link  of  the  maximal  spanning  tree  satisfies: 

(3)  *ip : Mln  (fij’  . V 

* 

f^  will  denote  the  maximal  flow  of  the  link  and  is  equal  to  its  capacity. 
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i 

where  f^,  f^,  ....  fQp  fora  the  unique  path  (X^,  ...»  Xop) 

connecting  nodes  jL  and  j>  within  the  tree.  For  If  the  Inequality  did  not 
hold,  the  smallest  link  of  the  path  could  be  removed  and  the  direct  link 
flp  substituted  to  form  a  tree  with  larger  value. 

From  (2)  and  (3),  we  conclude: 

(4)  f.  -  Min  <f..,  ...»  f  )  . 

ip  ij  op 

As  there  arc  only  n-1  links  in  the  tree,  we  have  shown  that  there  are 
at  most  n-1  different  values  of  the  maximal  flows. 

In  order  to  construct  the  maximal  spanning  tree,  we  introduce  the 
notion  of  a  condensed  network. 

Suppose  that  a  maximal  flow  problem  f>t  has  been  solved  between  the 
source  s,  and  the  sink  jt  by  locating  a  minimal  cut  (X,  X)  with  s.  in 

X  and  t.  in  X  .  To  find  f ^  where  nodes  JL  and  J.  are  at  the  same 

side  of  the  cut  (X,  X)  ,  say  in  X  ,  we  create  a  condensed  network  Q  by 

shrinking  all  the  nodes  of  X  into  a  single  node  to  which  all  the  arcs  of 

the  minimal  cut  are  attached.  Several  arcs  joining  the  same  pair  of  nodes 
can  be  replaced  by  a  single  arc  whose  capacity  is  given  by  the  sum  of  the 
capacities  of  the  arcs  composing  it. 

Lemma: 

The  maximal  flow  between  two  ordinary  nodes  ,1  and  in  the 

condensed  network  Q  is  numerically  equal  to  the  flow  f ^  in  C  . 

Proof : 

in  G  is  never  greater  than  the  flow  between  i.  and  J  in  Q 
because  shrinking  the  nodes  of  X  in  G  to  a  single  node  X  in  Q  is 


- l-W 


6 


equivalent  to  making  all  arcs  joining  the  nodes  of  X  of  infinite  capacity. 

On  the  other  hand  f ^  in  G  is  never  smaller  than  the  flow  in  Q 
because  if  the  f ^  in  Q  has  any  flow .  y  from  a  node  in  X  to  X  , 
there  is  a  set  of  arcs  from  X  to  X  in  the  original  network  with  a  total 
capacity  of  at  least  y  •  This  statement  can  be  made  about  every  arc 
connecting  the  set  X  with  the  node  X  . 

The  Gomory  and  Hu  procedure  for  constructing  the  tree  is  the  following: 

Take  any  two  nodes  and  find  a  minimal  cut  (X,  X)  separating  then.  We 
represent  this  by  two  nodes 


f 


1 


connected  by  a  link  with  the  value  of  the  cut.  In  X  are  listed  all  the 
nodes  of  the  set  X  and  in  X  the  others.  We  now  repeat  the  process. 

Choose  two  nodes  in  X  (or  two  in  X)  and  solve  the  flow  problem  in  the 
condensed  network  in  which  X  (or  X)  is  a  single  node.  The  resulting 
cut  has  a  value  f ^  and  is  represented  by  a  link  connecting  the  two  parts 
into  which  X  is  divided  by  the  cut,  say,  X^  ,  and  X2  .  X  is  attached 
to  X^  ll;  it  is  in  the  same  part  of  the  cut  as  X^  ,  to  X2  if  it  is  in 
the  same  part  as  X2  • 

The  procedure  is  continued  in  the  same  way  until  each  condensed  node 
consists  of  only  one  node. 

Examples  of  this  procedure  are  given  in  subsequent  illustrations  of  the 
algorithms  of  this  paper,  which  use  the  Gomory  and  Hu  method  as  a  subroutine. 


Lemma; 


The  maximal  flow  between  any  two  nodea  of  the  original  network  la  a  imply 
the  min  (f^,  f^*  •••»  the  •*rLes  of  links  of  the  tree  connecting 

the  two  nodes. 

We  now  prove  in  an  indirect  way  that  the  tree  constructed  as  above 
corresponds  to  the  maximal  spanning  tree  of  the  modified  network. 

We  have  f^  <  min  (f^,  fj . fr>  ;  for  each  f^  on  the  path 

connecting  nodes  jL  and  J.  ia  attached  to  a  cut  separating  JL  and  J,  . 

The  following  figures  illustrate  odr  discussion. 


numbers  in  arcs  are  capacities  and 
curve  lines  are  the  cuts 


corresponding  final  tree 


The  proof  of  the  reverse  Inequality  is  more  difficult.  It  will  be 

shown  that,  at  any  stage  of  the  construction,  if  a  link  of  capacity  J[ 

joins  nodes  X  and  7  in  the  tree,  then  there  is  an  x  in  X  and  a  y 

in  T  such  that  f  *  f  . 

xy  - 

The  proof  la  by  induction.  It  is  clearly  true  at  the  first  stage. 

Let's  suppose  it  is  true  in  the  rC**  stage  and  we  will  see  that  is  still  true 
in  the  next  one.  Consider  X  and  Y  ,  two  contiguous  condensed  nodes  of 
the  tree  at  the  r**1  stage  (t  links,  r+1  condensed  nodes)  and  xcX  and 
ycY  and  f  ■  f.  ■  capacity  of  the  link  (X,  Y)  .  Suppose  Y  has  more 

*y 

than  1  node,  at  least  has  nodes  s,  and  t.  ,  where  the  possibility  than  any  . 
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one  of  this  being  equal  to  £  is  not  excluded.  In  the  next  stage  the  set 

Y  divides  In  Y^  end  Yj  with  £  in  Y^  and  t_  In  Y£  .  We  may 

assume  that  x  Is  attached  to  Y^  . 

Now  s  and  t  are  the  two  nodes  such  that  f  _  -  f'  is  the 
—  —  st 

capacity  of  the  link  •  With  respect  to  t_  ,  there  are  two  easca. 


If  £  is  in  Y^,  f  is  still  equal  to 

show  that  f  ■  f  . 

xs  — 


By  Lemma  1,  condensing  Yj  to  a  single 
not  affect  f^  .  Denoting  by  bars  the 
we  have: 


_f  .  If  ^  is  in  Yj  ,  we  shall 


node  in  the  original  network  does 
maximal  flows  in  the  condensed  network. 


f 

f 

f 

f 


m 

f 

xs 

xs 

> 

f 

xy 

xy 

■ 

CD 

yt 

> 

f  * 

ts 

- 

ts 

f 


f'  . 


By  Theorem  2,  f  >  Min  [f  ,  f  .  ,  f _  ]  and  hence  f 
7  ’  xs  xy  yt  ts  xs 

But  f*  .>  f  because  a  cut  of  capacity  f'  separates  3C 
f  ■  f  ,  as  was  to  be  shown. 

XS  — 


f  >  Min  [f  ,r  ]  . 
xs  - 

Y  .  Thus 


and 


In  consequence,  the  capacities  of  the  links  in  the  final  tree  actually 


represent  maximal  flows  between  its  adjacent  nodes. 
Now  from  formulas  2  and  3  of  Theorem  2,  we  get: 


fip  -  Mln  (fij»  fjk*  fop) 

flp  5  Min  (fij*  fjk*  •••»  fop> 


where  i  j ,  Jk,  op  ,  are  the  links  of  the  tree  joining  i.  and  j>  , 

and  we  get  f.  -  Mln  (f,.,  ....  f  )  showing  that  the  final  tree  obtained 
ip  lj  op 

above  corresponds  to  the  maximal  spanning  tree  defined  earlier.  Because 
of  the  way  it  is  obtained,  we  will  call  this  tree  the  cut-tree. 

In  the  next  section,  we  discuss  the  Influence  on  the  cut-trae  of  a 
change  in  the  capacity  of  1  arc. 
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SECTION  II 

Sensitivity  Analysis  of  Multltermlnal  Flow  Networks  Under  Variation  of  the 
Capacity  of  1  Arc 

In  this  section  an  algorithm  is  developed  for  the  effect  on  the 
aaxlmal  flow  between  any  2  nodes  of  the  reduction  or  Increase  in  the  capacity 
of  any  arc.  The  capacity  of  all  arcs  except  arc  £  (where  arc  a  Is  any  of 
the  arcs  of  the  network)  will  be  considered  constant ,  and  arc  £  will  have 
a  capacity  C(a)  -  £  where  c  Is  a  parameter  that  can  take  all  the  values 
in  the  range  (0  ;  C(a)  ]  . 

Our  goal  is  to  find  the  relation  between  the  value  of  £  and  the 
resulting  cut-trees  for  an  undirected  network. 

There  are  ranges  of  the  value  of  t_  ,  let  us  say  for  all  £'s  between 
Cj  and  e 2  ,  for  which  one  cut-tree  is  enough  to  describe  the  maximal  flows 
between  all  pairs  of  nodes;  moreover,  there  is  a  finite  number  of  cut-trees 
that  successively  correspond  to  the  variation  of  the  value  of  £  from  0  to 
C(a)  . 


h— 

i 

- 1 — 1 - 

0 

i 

i 

< 

• 

!  C(a) 

1 

• 

« 

1 - 

« 

$ 

1 

- 1 - 

1 

• 

« 

- 1 - 

Cut-tree  Cut-tree  Cut-tree 
1  2  3 


Algorithm 

Step  1: 

Construct  the  cut-tree  of  the  network: 

a)  In  each  step  of  this  construction  put  an  initial  capacity  of  zero 
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to  the  arc  a.  .  Determine  the  maximal  flow  between  the  two  nodes 
currently  designated  source  and  sink, 
b)  Reset  the  capacity  of  arc  11  to  C(a)  -  c  and  continue  maximizing 
the  flow  until  it  results  in  a  cut-set.  The  value  of  the  flow 
f(X,X)  through  each  cut-set  (X,X)  can  be  an  expression  of  one  of 
these  three  types: 

1)  A  constant  a  (this  means  that  the  flow  didn't  change  from 
a)  to  b)  and  the  cut-set  is  not  through  arc  a)  . 

II)  A  function  of  e  as:  y-e  where  y  is  a  constant.  (Cut¬ 
set  through  arc  a  .) 

Ill)  A  constant  6  if  c  <  k  .  (The  maximal  flow  did  change  from 
a)  to  b)  and  the  cut-set  is  not  through  arc  a_  for  c  <  k  .) 

We  record  the  value  of  k  associated  with  each  6  . 

Step  2: 

a)  If  the  cut-tree  doesn't  have  any  k  it  will  represent  the  maximal 

*  a 

flows  between  any  pair  of  nodes,  for  0  <  e  <  C(a  )  where  C  (a) 
will  represent  the  current  value  of  C(a)  ,  the  value  of  the 
capacity  in  the  specific  cycle.  (Initially  C  (a)  ■  C(a).) 

b)  If  there  are  some  values  of  k  ,  take  the  smallest  of  them  and 
recompute  the  maximal-flow  between  the  two  nodes  of  the  link.  Two 
things  can  happen:  Either  the  value  of  such  k  remains  the  same 

in  which  case  follow  Step  (3)  or  it  is  Increased  in  which  case  it  has 
to  be  compared  to  the  others  and  if  it  is  still  the  smallest,  go  to 
Step  (3),  otherwise  go  back  to  Step  (2),  phase  b. 

Step  3: 

* 

After  all  the  computations  in  Step  (2) ,  we  let  k  -  smallest 

(k  recomputed}. 

*  * 

If  k  >  C  (a),  this  is  the  last  cycle  and  as  in  Step  (2),  the  tree 
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represents  the  maximal  flows  for  all  the  positive  values  of  the 
current  capacity  of  arc  a  ;  otherwise  replace  the  capacity  of  arc 

a  ■  (C  (a)  -  c)  by  ((C  (a)  -  k  )  -»  e)  where  (C  (a)  -  k  )  will  be 

* 

the  current  value  *f  C  (a)  for  the  next  cycle.  Go  to  Step  (1). 

Note : 

There  is  an  example  application  of  the  algorithm  at  the  end  of  the 
paper. 

The  network  used  is  the  same  one  used  by  Elmaghraby  to  facilitate  a 
comparison  of  our  method  and  his. 

Justification  of  the  Algorithm 

In  the  construction  of  the  cut-tree,  the  maximum  flow  between  two  nodes 
1  and  J  is  determined  In  one  of  two  ways:  Either  directly  by  a  standard 
computation  of  the  maximum  flow  between  these  nodes  or  indirectly  as  a 
result  of  computing  maximal  flows  between  other  nodes.  If  the  maximum  flow 
Is  found  by  a  direct  computation,  sending  the  maximum  possible  amount  of 
flow  through  arcs  other  than  arc  a  ,  and  only  subsequently  using  the  arc  a  , 
we  Indeed  get  the  maximum  value  of  k  for  which  the  expression  of  the  link 
represents  the  maximum  flow  between  the  two  nodes. 

If  the  maximum  flow  between  nodes  i^  and  Is  determined  by  computating 
of  the  flow  between  nodes  r_  and  £  where  jr  of  a_  or  both  are  different 
from  1_  and  ,  and  if  the  capacity  of  the  link  (1,J)  is  of  the  form 
B(k)  ,  this  value  of  k  could  have  been  constrained  by  any  other  link  of 
the  type  II,  and  so  a  function  of  e  .  In  this  case  it  will  be  detected  in 
Step  (2)  of  the  algorithm. 

An  example  will  help  to  understand  this  situation.  The  numbers  in  the 
trees  and  links  are  the  respective  capacities. 
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Let's  suppose  that  after  a  first  step  in  developing  a  cut-tree,  we 
arrive  at  this  situation: 


reduced  network  cut-tree 


If  we  choose  nodes  (b)  and  (?)  for  the  next  step,  we  get:  (The  numbers 
in  heavy  black  are  the  flows). 


Maximum  flow  =  6  +  8  +  (10  -e)  *  24  -  c  . 
Now  solving  for  nodes  and  (?) 
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Maximum  flow: 


15+6+2  [8]  -  23  (8] 

(That  is,  maximum  flow  is  23  if  e  <  8  .) 

There  is  only  one  link ,  namely  the  one  connecting  nodes  and  Q)  that 
has  an  associated  k  value:  k  *  8  .  This  link  was  computed  directly 
because  we  did  solve  the  maximum  flow  between  and  . 

Now  we  will  see  what  happers  if  that  link  is  obtained  in  another  way. 
Beginning  again,  we  now  look  for  the  maximum  flow  between  ^  and  . 


Maximum  flow:  8+6+9  [1]  *  23  [1] 

Between  @  and  ®  we  get  the  same  maximum  flow  as  before.  So  the  cut-tree 
will  now  be: 
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We  see  that  in  this  case  we  have  k  =  1  ,  but  the  link  (?)  -  was  not 

directly  computed,  (we  computed  links  (b)  -  and  -  (?).)  So  we  can 

compute  ©  -  ©  to  check  the  value  k  •  1  and  (as  before)  we  get  k  -  8  . 

Applications 

Following  are  some  problems  to  which  our  algorithm  can  be  applied. 

1.  Electricity  line  network  of  a  country  under  partial  failure  or  reparation 
of  a  line  between  two  cities. 

2.  Similar  problems  with  commodities  like  water  or  gas  in  a  city,  or  flows 
of  water  for  the  event  of  fire  in  a  "dry"  city. 

3.  City  of  one-way  (or  almost  one  way)  freeways  where  the  direction  of  the 
freeways  can  be  reversed,  and  where  the  people  go  to  different  places  at 
specific  times  and  the  optimal  number  of  lanes  of  a  new  highway  is  under 
study . 

4.  Constructing  roads  in  a  war  problem  "Jungle": 

Given  a  set  of  detachments  (nodes)  defending  an  area  and  very  limited 
means  of  transportation  (capacities),  find  a  maximum  fl'"-*  of  soldiers 
between  the  detachment  that  happens  to  be  idle  to  the  one  under  attack. 


Cut  set  between 


and  Q>). 


(reduced  network) 


cut-tree. 


Cut  set  between  and  . 
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Cut  set  between  and  . 
(reduced  network) 


(27-e)  does  not  Interfere 

because  maximum  (c)  **  13 


<? 


The  link  (3,6)  with  value  37(3]  is  the  one  that  has  the  minimum  k  ,  k  -  3  . 
As  link  (3,6)  ,  was  computed  directly,  it  is  not  possible  to  improve  that 
value  and  still  have  the  same  cut  tree.  This  cut  tree  represents  the  maximum 
flows  of  the  original  network  for  0  <  e  <  3  . 

Now,  we  replace  the  capacity  of  arc  (2,5)  by  13-3-e  »  10-c  and  begin 


again. 


Cut  set  between  (?)  and  (7). 


(reduced  network) 


Cut  tree. 


Cut  set  between  ^)  and 


Cut  set  between  (2)  and 


Cut  tree. 


0  0 


23 


24-e 

O19  33_C  A  37-e  A  n[l]r\ 

- ©— — © - ©— ^ 


The  only  link  with  k  >  0  is  link  (6,7)  with  k  *  1  ;  as  link  (6,7)  is  not 
the  result  of  the  computation  of  the  cut  set  between  ®  and  but  instead 
it  was  of  (j)  -  (7)  and  ®  ,  we  have  to  compute  the  maximum  flow  between 


Cut  tree. 


24 


Third  cycle 

Original  network  with  the  capacity  of  arc  (2,5)  reduced  to  2-r.  . 


Cut  set  between  and  . 


Cut  tree. 


29-e 


Cut  set  between  Q)  and  (?) 


Cut  tree. 

C 


Cut  set  between  and 


Cut  tree. 


G> 


19 


23 -c 


<Z> 


16- 


We  did  not  get  any  value  of  k  ,  so,  this  cut  tree  represents  the  original 
networx  for  the  maximal  flows  between  any  two  nodes  for  11  <  e  <  13  .  i.e., 

e  *  11  In  the  original  network  corresponds  to  e  -  0  in  this  cut  tree;  e  ■ 

in  the  original  network  corresponds  to  e  ■  2  in  this  cut  tree. 

i.e. ,  maximal  flow  between  and  Q)  is  23  regardless  of  c  , 

"  "  ©  "  "  ©  "  the  minimum  of 

(23  and  25-c)  . 

"  "  is  the  minimum  of 

(23-e  ,  29 -e ,  25-c  and  19)  , 


etc. 
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SECTION  III 

SENSITIVITY  ANAI/'Sl S  OF  MULTITERM  1 NAh  FLOW  NETWORKS  UNDER  VAK1 AT  I ON  OF  THE 
CAPACITY  OF  SKVEKAI."  ARCS 

We  generalized  the  results  of  the  previous  section  by  allowing  several 
of  the  arcs  of  the  network  to  be  linear  functions  of  a  parameter.  Further 
generalizations  (several  parameters  or  nonlinear  functions  of  them)  will 
not  be  attempted  in  the  present  work. 

The  algorithm  is  very  similar  to  the  one  discussed  in  the  first  part. 
Now  we  have  arcs  with  capacities  of  three  types: 


+  a 


ij 


x 


>  0  . 


We  will  try  to  find  the  sequence  of  cut-trees  that  represent  the  maximal 
flows  between  any  two  arcs,  when  x  varies  over  the  real  numbers  between  0 
and  some  arbitrary  maximum  x  . 

When  the  capacity  of  an  arc  becomes  0  ,  we  simply  erase  such  an  arc  not 
having  defined  negative  capacities. 


The  Algorithm 


Step  1: 

Following  the  Gomory-Hu  procedure,  we  construct  the  cut-tree  of  the 
network,  in  each  step,  using  the  labeling  algorithm  of  Ford-Fulkerson  to 
maximize  the  flow  between  nodes  r  and  s  (generic  nodes)  . 

As  a  result  of  the  above,  we  will  have  in  each  arc  (i,j)  a  flow* 
f^  and  a  residual  capacity  r^  representing  the  difference  between  the 
capacity  and  the  flow  through  the  arc.  Also,  we  will  have  a  cut-set  and  the 
value  of  its  capacity. 

f^j  can  be  of  one  of  three  types:  f^  ;  f  +  °ijx  anc*  “  aijx 
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<a  •  “u  i 0  • 

can  be  of  one  of  three  types:  +  Bjj  ;  -  B^x  • 

ru  •  6u  : 0  • 

For  the  flows  and  residual  capacities  of  the  types  f  -  a^x  and 
r^j  -  8jjx  •  ve  compute  the  maximum  value  of  x,  call  It  k  ,  that  we  can 
have  before  any  flow  or  capacity  becomes  negative. 

Three  cases  can  occur: 

1)  There  Is  no  value  of  k  ,  i.e.,  no  flows  go  through  arcs  that  can 
be  affected  by  x  .  In  this  case,  compute  the  next  link. 

2)  The  value  of  k  Is  greater  than  the  maximum  x  ,  we  want  to 
consider.  Go  to  the  next  link. 

3)  The  value  of  k  is  greater  than  0  but  smaller  than  maximum  x  . 
In  this  case,  reset  the  value  of  the  capacities  that  depend  on 

x  tc  :  +  (k  +  a^)x  arid  -  (k  +  a^j)x  • 


Compute  again  the  same  maximum  flow  in  this  network  with 
modified  capacities. 

If  both  the  cut-set  and  the  expression  of  maximum  flow  don't  change, 
we  get  a  new  value  of  k  that  should  be  added  to  the  previous  k  , 
and  if  this  value  is  still  smaller  than  maximum  x  we  repeat  again 
and  again  +  until  either 


If  the  network  is  not  too  complicated  an  intelligent  assignment  of  flow3 
will  frequently  yield  the  optimal  values  in  the  first  try  without  requiring 
successive  adjustments.  In  the  example  that  follows  only  two  or  three 
repetitions  were  required  at  each  stage.  The  value  of  the  residual 
capacity  is  written  beside  the  flow  only  when  necessary  for  understanding. 
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we  are  in  Case  2.  or  either  the  cut-set  or  the  expression  of  the 
maximum  flow  change  in  which  case  we  have  finished  the  computation 
of  the  link. 


Step  IT : 

Applies  when  we  have  completed  the  computation  of  the  n-1  links. 

Now  we  have  the  cut-tree. 

Take  the  link  with  the  smallest  value  of  k  .  If  this  value  is  larger 
than  maximum  x  ,  we  have  finished,  otherwise,  if  the  link  was  computed 
directly  go  to  Step  3;  if  it  was  not  computed  directly  compute  again  the 
value  of  this  link  getting  the  same  value  of  k  or  a  larger  one.  Go  back  to 
Step  2. 

Step  III: 

The  value  of  k  obtained  in  Step  2  gives  the  range  for  what  the  last 
tree  is  valid. 

I.e.,  the  first  tree  is  valid  for  0  <  x  <  k  ,  the  second  for 
<  x  <  (k^  +  kj)  etc. 

Reset  the  values  of  the  capacities  that  are  functions  of  x  to 


f  Ctj  +  ka4J  + 

cij  +  “ijx 

and 

Instead  of 

cij  -  koij  ■  “ijx 

'  * 

Clj  •  °1JX 

go  back  to  Step  I  for  a  new  cycle. 

Justification  of  the  Algorithm 

In  Step  I,  we  maximize  the  flow  between  two  nodes.  Thi9  can  be 
accomplished  in  many  ways.  Each  way  leaves  some  residual  capacities  and 


flows  as  a  function  of  x  that  tend  to  zero  with  the  increase  of  the  value 
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I 

of  x  ,  yielding  an  upper  bound  k  on  the  value  of  x  . 
Changing  the  capacities  of  the  arcs  of  the  types 


C1J +  °ljx 

to 

(cij  +  kaij)  +  aijx 

CJJ  -  °1JX 

(cij  *  kaij>  -  aijx 

we  can  change  the  amount  of  flow  carried  by  the  different  arcs  to  get  the 
minimum  restriction  and  hence  the  maximal  value  of  k  . 

In  Step  II,  as  in  the  previous  algorithm  (and  for  the  same  reasons), 
we  compute  again  the  k  corresponding  to  the  link  with  minimum  k  ,  if  it 
was  not  computed  directly. 

Important 

In  any  direct  computation  the  value  of  k  can  get  larger  or  stay  the 
same,  but  never  get  smaller,  because  it  is  limited  by  all  the  restrictions 
imposed  by  the  cut-sets  that  were  computed  indirectly  to  get  the  link  we 
are  considering. 

When  we  compute  the  link  directly,  fewer  restrictions  result  and  the 
maximum  value  of  a  function  cannot  diminish. 

Examples  of  Applications 

1)  The  network  consists  of  one  way  highways,  some  of  asphalt  and 
some  of  concrete. 

The  capacity  of  the  highways  will  be  given  the  number  of  lanes, 
and  the  speed  of  the  cars  on  the  asphalt  lanes  decreases 
proportionally  (for  some  range)  to  the  amount  of  rain.  Thus,  the 
capacity  of  a  four  lane  concrete  highway  is  4,  while  the  capacity 


of  a  4  lane  asphalt  highway  is  4  -  4x  ,  etc. 

2)  The  network  consists  of  electric  lines  going  from  power  plaits  to 
the  cities.  Some  of  the  plants  are  hydraulic  plants  and  others  an 
thermal.  Because  of  increase  of  coal  costs  it  is  under  const  derat  •  i 
to  decrease  the  capacity  of  the  coal  plants  and  to  increase  the 
capacity  of  the  hydraulic  ones.  (In  this  case  artificial  arcs  with 
the  capacities  of  the  plants  have  to  be  added  to  the  network.) 

Comments  for  the  Second  Algorithm 

The  problem  solved  by  the  second  algorithm  was  taken  from  an  unpublishe< 
paper  by  Elmaghraby.  The  example  presented  is  also  his,  to  facilitate  an 
eventual  comparison  of  methods. 

Step  I  of  the  present  algorithm  is  essentially  similar  to  the  one  in 
his  paper. 

Step  II,  which  differs  from  the  approach  of  Elmaghraby,  provides  a  way 
to  bypass  the  construction  of  superfluous  cut-trees. 
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Example 

Or ip; Inal  network 

Sensitivity  analysis  will  be  done  for  X  in  the  range  for  this  m  1 1 i  — 
terminal  network.  0  <  X  <  6.66 


Construction  of  the  first  cut-tree  or  (first  cycle). 


changes. 


No  iteration  was  done  because  it  is  evident  that  for  X  >  7 
t 


csx 

for 

cs2 

for 


the  cut-set 

Flow  »  19 
is  the  cut-set 
X  <  7  . 
is  the  cut-set 
X  >  7  . 


The  same  thing  happens  most  of  the  time  throughout  this  example. 
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MV] 


k.4+^4.25 


0 


f;,3  .wJrWJSLQ 
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Link  (2,6)  has  the  smallest  value  of  k  ,  k  -  1  .  As  this  link  was 
not  computed  directly  but  instead  of  it  were  (2,3)  and  (3,6),  we  have  to 
check  its  value. 


rs 


As  we  see,  the  value  of  k  ,  has  been  improved  to  k  •  2  .  It  happens 
that  it  is  still  the  smallest,  and  this  time  it  has  been  computed  directly, 
so  we  cannot  improve  any  more. 

This  cut-tree  (the  last  one)  represents  the  maximum  flows  in  the 
network  for  values  of  x  :  0  <  X  <  2  .  The  next  cut-tree  will  represent 

the  maximum  flows  for  values  of  X  larger  than  2  .  The  upper  limit  will 
be  found  in  the  same  way  we  got  the  upper  limit  2  for  the  first  cut-tree. 


first  iteration 


5  -  X  -  0 


second  iteration 


***121  /O 


Flow  »  27 


necessary  to  compute  it  again. 


it  is  definitive. 


A3 


7.5-X-O  - >  x-7.5 


23. 5-3X[3] 


■G> 


28.5-X[7.5] 


o 


2711.51 


<D 


No  restriction 


19 


;.M-;3,5~3xm-{^)-28'5~x|7-sl  0 


Flow  -  19 
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SECTION  IV 


COMMENTS,  COMPAR1 SONS 

Professor  Salah  E.  Elmaghraby  of  the  Department  of  Iiuhistrial 
Administration  of  Yale  University  has  two  algorithms  [2]  and  (4) 
that  solve  the  same  problems. 
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